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ABSTRACT 

For inviscid, rotational accretion flows, both isothermal and polytropic, a simple dy- 
namical systems analysis of the critical points has given a very accurate mathematical 
scheme to understand the nature of these points, for any pseudo-potential by which 
the flow may be driven on to a Schwarzschild black hole. This allows for a complete 
classification of the critical points for a wide range of flow parameters, and shows that 
the only possible critical points for this kind of flow are saddle points and centre-type 
points. A restrictive upper bound on the angular momentum of critical solutions has 
been established. A time-dependent perturbative study reveals that the form of the 
perturbation equation, for both isothermal and polytropic flows, is invariant under the 
choice of any particular pseudo-potential. Under generically true outer boundary con- 
ditions, the inviscid flow has been shown to be stable under an adiabatic and radially 
propagating perturbation. The perturbation equation has also served the dual pur- 
pose of enabling an understanding of the acoustic geometry for inviscid and rotational 
flows. 
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INTRODUCTION 



Multi-transonicity in axisymmetric accretion becomes possible for a certain range of values of the specific angular momentum 
of the flow, when the number of critical points in the accreting system exceeds one, unlike in the case of spherically symmetric 
accretion. A substantial body of work o ver the past many years, has argued well the case for multi-transonicity in axisymmetric 
accretion lAbramowicz fc Zureklll98lt [Fukuelll987l: IChakrabartJll99(]|: IXafatos fc YandllQQjJ: IVaiig fc Kafa,tojll995l: IPariev 



h99g': 'Lasota fc Abramowic j |l997l: ILu et alJll997l: IPeitz fc Appllll997l: lOaslbooi iBarai et ahlbood : lOajboO^ lAbraham et al 
12006 : Das ct al. 2006). Transonicity would imply that the bulk velocity of the flow would be matched by the speed of acoustic 
propagation in the accreting fluid. In this situation a subsonic-to-supersonic transition or vice versa can take place in the 
flow either continuously or discontinuously. In the former case the flow is smooth and regular through a critical point for 
transonicity (more specially, this ca n be a sonic point, a t which the bulk flow exactly matches the speed of sound), while in the 
latter case, there will arise a shock JChakrabartilll98gl) . The possibility of both kinds of transition happening in an accreting 
system is very much real, and much effort has been made so far in studying these phenomena. For accretion on to a black 
hole especially, the argument that the inner bou ndary condition at the event horizon will lead to the exhibition of transonic 
properties in the fiow, has been established well l)Chakrabartilll99(l l . However, for all the importance of transonic flows, there 
exists as yet no general mathematical prescription that allows for a direct understanding of the nature of the critical points of 
the flow (and the physical solutions which pass through them) without having to take recourse to the conventional approach 
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of numerically integrating the governing non-linear flow equations, and then looking for any physically meaningful transonic 
behaviour under a given set of flow parameters. 

This paper purports to addr ess that particular issu e. To do so it has been necessary to adopt mathematical methods from 
the study of dynamical systems Jjordan fc Smithlll99gl) . In accretion liter ature one would come across some relatively recent 
work s which had indeed made use of the techniques of dynamical systems llRav fc BhattacharieelboO^ : lAfshordi fc Paczvhskil 
l2003l) . but the present treatment, it is being claimed, is much broader and more comprehensive in its scope and objectives, 
than those reported previously. Through this approach a complete and mathematically rigorous prescription has been made for 
the exclusive nature of the critical points in axisymmetric pseudo-Schwarzschild accretion, and, in consequence, the possible 
behaviour of the flow solutions passing through those points and their immediate neighbourhood. In conjunction with a 
knowledge of the boundary conditions for physically feasible inflow solutions, this makes it possible to form an immediate 
qualitative notion of the solution topologies. What is more, this study has been shown to set a more restrictive condition on 
sub-Keplerian flows passing through the critical points. 

Beyond this stage, time-dependence has been brought into the flow equations. This has been easy because, although 
this study is devoted to black hole accretion, the Newtonian construct of space and time has been preserved, with the 
general relativistic effects of strong gravity in the vicinity of a Schwarzschild black hole having be en represented through a 
modified pseudo-Newtonian potential, for whic h, various prescriptions exist in accretion literature JPaczvhski fc Wiitalll98ol : 
iNowak fc Wagoneilll99ll lArtemova etai]|l996l) . Thus it has been possible in this Newtonian framework, to carry out a time- 
dependent stability analysis of the physically well-behaved and relevant stationary solutions in a straightforward manner. The 
stability of these solutions has been shown to be dependent on the asymptotic conditions on large length scales of the accretion 
disc, and since these conditions are the same for the choice of any pseudo-Newtonian potential, it has been argued that the 
flow solutions will remain stable under linearised perturbations in all cases. Another interesting aspect of this perturbative 
study in real time has been that the equation of motion for the propagation of the perturbation yields a metr ic that is entirely 
identical to the metric of an analog acoustic black hole obtained through very different means llVisseilliggj) . 

One important feature of the mathematical treatment in this paper is worth stressing. The entire study — on the nature 
of the critical points as well as on the time-dependent stability of the ffow — has been carried out in terms of a general 
potential only. Therefore, the mathematical results will hold equally well under the choice of any particular potential that 
will drive the flow — both isothermal as well as a general polytropic — on to the black hole. 

Finally, the range of applicability of the mathematical methods of dynamical systems goes much beyond a gravitational 
system deflned purely within the Newtonian framework. In a rigorously general relativistic stationary flow, described variously 
by either the Schwarzschild metric or the Kerr metric, it is eminently possible to carry out a mathematical analysis in the 
same principle in which it has been presented in this work. These general relativistic cases will be reported separately, while 
the present work may be considered to be the flrst in a series in which the usefulness of a dynamical systems approach in 
astrophysical problems has been considered at length. 



2 THE EQUATIONS OF THE FLOW AND ITS FIXED POINTS 

It is a standard practice to consider a thin, rotating, axisymmetric, inviscid steady flow, with the condition of hydrostatic 
equilibrium imposed along the transverse direction ijMatsumoto et al]ll984l : IPrank et ahlbnO^) . The two equations which 
determine the drift in the radial direction are Euler's equation, 

dr p ar r-^ 
and the equation of continuity, 

A i^p^rH) = (2) 

in which, 0(r) is the generalised pseudo- Newtonian potential driving the flow (with the prime denoting a spatial derivative), 
A is the conserved angular momentum of the flow, P is the pressure of the flowing gas and H = H{r) is the local thickness of 
the disc JPrank et al.ll2002ll . respectively. 

The pressure, P, is prescribed by an equation of state for the flow llChandrasekhailll939l) . As a general polytropic it is 
given as P = Kp^ , while for an isothermal flow the pressure is given by P = pnT/ fimu, in all of which, K is & measure of the 
entropy in the flow, 7 is the polytropic exponent, k is Boltzmann's constant, T is the constant temperature, mn is the mass 
of a hydrogen atom and fi i s the reduced mas s, respectively. The function H in Eq.ljJJ will be determined according to the 
way P has been prescribed ijPrank et alJl20o3l . while transonicity in the flow will be measured by scaling the bulk velocity 
of the flow with the help of the local speed of sound, given as Cg = {dP / dpY^^ . In what follows, the flow properties will be 
taken up separately for the two cases, i.e. polytropic and isothermal. 
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2.1 Poly tropic flows 



With tiie polytropic relation specified for P, it is a straightforward exercise to set down in terms of the speed of sound, Cs, a 
first integral of Eq.Q as, 



— + nCs + > 



ir) + 



2r2 



(3) 



in which n = (7 — and the integration constant £ is the Bernoulli constant. Before moving on to find the first integral 
of Eq.||2l it should be important to derive the functional form of H. Assumption of hydrostatic equilibrium in the vertical 
direction deliver this form to be 

1/2 



H : 



with the help of which, the first integral of Eq.l|5J could be recast as 



2(2„+i) y^'r^ 



(4) 



(5) 



where M = (7,/^ )"m Jchakrabarttil99d) with m, an integration constant itself, being physically the matter flow rate. 

To obtain the critical points of the flow, it should be necessary first to differentiate both Eqs.Q and JSJ, and then, on 
combining the two resultant expressions, to arrive at 



(2 ^2 2\ d , 2n 



2v' 



(6) 



with 0^ = 2(7 + The critical points of the flow will be given by the condition that the entire right hand side of Eq.l|HJ 

will vanish along with the coefficient of d{v^)/Ar. Explicitly written down, following some rearrangement of terms, this will 
give the two critical point conditions as, 

-1 



fl2 2 
P ^sc 



rc(p'{rc) ■ 



(7) 



with the subscript c labelling critical point values. 

To flx the critical point coordinates, Vc and rc, in terms of the system constants, one would have to make use of the 
conditions given by Eqs.JTJ along with Eq.J^I, to obtain 

-1 



27 



7- 



1 



)2 

rc(j}'irc) - -j 



(8) 



from which it is easy to see that solutions of rc may be obtained in terms of A and £ only, i.e. rc = fi{X,£). Alternatively, Vc 
could be fixed in terms of A and A4. By making use of the critical point conditions in Eq.J^ one could write 

, 2(n + l) 



A^' 

rc0'(rc) 2 



<t>'(rc) 



(9) 



with the obvious implication being that the dependence of rc will be given as Tc = f2{\,M). Comparing these two alternative 
means of fixing rc, the next logical step w ould be to say that for the fixed points, and for the solutions passing through them, 
it should suffice to specify either £ or (|Chakra bartilll99nl) . 



2.2 Isothermal flows 

For isothermal flows, one has to go back to Eq.Q and use the linear dependence between P and p as the appropriate equation 
of state. On doing so, the first integral of Eq.Q is given as 

f- + c!\np+^{r) + ^=C (10) 

with C being a constant of integration. For flow solutions which specifically decay out to zero at very large distances, the 
constant C can be determined in terms of the "ambient conditions" as C = Cg Inpoo. The thickness of the disc can be shown 
to bear the dependence 




with the help of which, the first integral of the continuity equation is given as. 
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(12) 



As has been done for polytropic flows, both Egs. llOt and 11211 are to be differentiated and the resuhs combined to give, 



^ - r<j!> + -Cs 3 - r — 



from which the critical point conditions are easily identified as, 

<^"(rc)" 



rc<t>'{rc) 2" 



(13) 



(14) 



In this isothermal system, the speed of sound, Cs, is globally a constant, and so having arrived at the critical point conditions, 
it should be easy to see that Tc and Wc have already been fixed in terms of a global constant of the system. The speed of sound 
can further be written in terms of the temperature of the system as Cs = QT^^^, where O = {k/ ^mu)^^'^ , and, therefore, it 
should be entirely possible to give a functional dependence for rc, as Tc — /a (A, T). 



3 NATURE OF THE FIXED POINTS : A DYNAMICAL SYSTEMS STUDY 



The equations governing the flow in an accreting system are in general first-order non-linear differential equations. There is 
no standard prescription for a rigorous mathematical analysis of these equations. Therefore, for any understanding of the 
behaviour of the flow solutions, a numerical integration is in most cases the only recourse. On the other hand, an alternative 
approac h could be made to t his question, if the governing equations are set up to form a st andard first-order dynamical 
system (|.Tordan fc SmithlllQQtt l. This is a very usual practice in general fluid dynamical studies ijBohr et al.lll99,4 l. and short 
of carrying out any numerical integration, this approach allows for gaining physical insight into the behaviour of the flows 
to a surprising extent. As a first step towards this end, for the stationary polytropic fiow, as given by Ea. |^. it should be 
neces sary to parametrise this equation and set up a coupled autonomous first-order dynamical system as jjordan fc SmithI 



dr' 



2«" 



// I ^ a2 2 
r(j) + -P 



dr / 2 

d7 = 



(15) 



in which r is an arbitrary mathematical parameter. With respect to accretion stud ies in particular, this kind of parametrisation 
has been reported before jRav fc Bhattacharieel2'oo3lAfshordi fc Paczviiskilboosl) . but here the possibility of a more thorough 
mathematical description of the nature of the critical points is being explored. 

The critical points have themselves been fixed in terms of the flow constants. About these flxed point values, upon using 
a perturbation prescription of the kind = -f Sv'^ , ci = c^c + Sc^ and r = rc + Sr, one could derive a set of two autonomous 
first-order linear differential equations in the 5r — Sv^ plane, with 5c^ itself having to be first expressed in terms of 5r and 
Sv'^ , with the help of Eq.@ — the continuity equation — as 



Sc. 

72 



7- 



5v^ 
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[3 .'^"('^c)] 









7+1 

The resulting coupled set of linear equations in 5r and Sv^ will be given as 



(16) 



1 IX 2^ 
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rc dr 
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(17) 



A = rc 



Trying solutions of the kind Sv^ 
Sv'^ and 5r — as 

4rc0'(rc)Cs, 



■/''"('•c) _ 0"(rc) 

■J exp(S7r) and 5r ~ exp(f2r) in Egs. llTL will deliver the eigenvalues Q, — growth rates of 



(7+1)^ 
where ^Y.if) = r'^(l>'{r) 



(7 - 1) ^ - 27 (4 ^) + 27S 1 
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47 + (7 - 1) ^ 27B 1 



(18) 



Multi-transonic pseudo-Schwarzschild accretion 5 



For isothermal flows, starting from Ea. ll3ll . a similar expression for the related eigenvalues may likewise be derived. The 
algebra in this case is much simpler and it is an easy exercise to assure oneself that for isothermal flows one simply needs to set 
7 = 1 in Eq. llBII . to arrive at a corresponding relation for f2^. However, it should be incorrect to assume that in this kind of 
study, one could always treat isothermal flows simply as a special physical case of general polytopic flows. For polytropic flows 
the position of the fixed points, under a given form of ^(r), will be determined either by Eq.Q or by Eq.Q. For isothermal 
flows the fixed points are simply to be determined from the critical point conditions, as Eas. H14|l give them, since Cs is globally 
a constant in this case. 

Once the position of a critical point, Vc, has been ascertained, it is then a straightforward task to flnd the nature of that 
critical point by using rc in Eg. 1181 . Since it has been discussed in Section |21 that is a function of A and T for isothermal 
flows, and a function of A and £ (or A4) for polytropic flows, it effectively implies that Q,^ can, in principle, be rendered as 
a function of the flow parameters for either kind of flow. A generic conclusion that can be drawn about the critical points 
from the form of Q,^ in Eg. llBH . is that for a conserved pseudo-Schwarzschild axisymmetric flow driven by any potential, the 
only admissible critical points will be saddle points and centre-type points. For a saddle point, > 0, while for a centre- type 
point, Q,'^ < 0. Once the behaviour of all the physically relevant critical points has been understood in this way, a complete 
qualitative picture of the flow solutions passing through these points (if they are saddle points), or in the neighbourhood of 
these points (if they are centre-type points) , can be constructed, a long with an impression of the direction that these solutions 
can have in the phase portrait of the flow (|,Tordan fc Smithlll999l) . 

A further interesting point that can be appreciated from the derived form of Q,^ , is related to the admissible range 
of values for a sub-Keplerian flow passing through a saddle point. It is self-evident that for this kind of flow, (A/Ak)^ < 
1 ^Abramowicz fc Zxireklll98lh . However, a look at Ea. ljlB^ will reveal that a more restrictive upper bound on A/Ak can be 
imposed under the requirement that Q,^ > for a saddle point, and this restriction will naturally be applicable to solutions 
which pass through such a point. This is entirely a physical conclusion, and yet its establishing has been achieved through a 
mathematical parametrisation of a dynamical system. 



4 NUMERICAL RESULTS 

It has been argued with the help of Ea. ll8ll . that the spatial coordinate of the critical points, Tc (which in its turn depends 
on the flow parameters), will determine the eigenvalues and the corresponding nature of a given critical point. This will 
also necessitate an explicitly given functional form of (j}(r). A pseudo-Schwarzschild flow is driven by a pseudo-Newtonian 
potential, which describes general relativistic effects which are most important for accretion disc structures in the vicinity of 
a Schwarzschild black hole. Introduction of such potentials will allow for the investigation of complicated physical processes 
taking place in disc accretion in a semi-Newtonian framework, by avoiding the difficulties of a purely general relativistic 
treatment. Through this approach most of the features of space-time around a compact object are retained and some crucial 
properties of analogous relativistic solutions of disc structures can be reproduced with a high degree of accuracy. 

In the present treatment, four such pseudo-Newtonian potentials have been chosen to make a graphical representation 
of both the qualitative and the quantitative features of the critical points with respect to a given set of fiow parameters {£, 
7 and A for polytropic fiows, and T and A for isothermal fiows). In general a potential has been represented by 4'i(r), with 
{i = 1, 2, 3, 4}. Written explicitly, each of these potentials is given as 

1 



Mr) = 



2 (r - 1) 



Mr) - 



l-A + 12fl 
2r V2r 



IN 1/2 

Mr) = 



r 



Mr) = (19) 

in all of which, the length of the radial coordinate, r, has been scaled in units of Schwarzschild radius, defined as — 
2GMbh/c^ (with Mbh being the mass of the black hole, G the universal gravitational constant and c th e velocity of light in 
vacuu m) . At various points of time, eac h of these p otentials has been intro duced in accretion literature bv lPaczvnski fc Wiital 
Jl980l) for (Ai . iNowak fc" Waeoneil il99lh for (f>2 , and lArtemova et al.' figg^ for 03 and (jn^ re specti vely. A compara t ive ov erview 
of the physical properties of these potentials has been given bv lPas fc Sarkarl j200ll) . lDa3 liooi) and bas et all tooi) . 

Considering polytropic flows first. Fig. Q gives the plots of £ versus A parameter space for all the potentials cpiir). The 
region marked O gives the values of £ and A corresponding to the single outer critical point. Similarly I indicates the region 
of the parameter space that will yield the lone inner critical point. Multi-transonicity is to be obtained for the values of £ 
and A contained within the wedge-shaped region, which is further subdivided into two regions, A (for accretion) and W (for 
wind). The area of multi-transonicity shifts according to the choice of a particular potential. 
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Figure 2. Variation of Q? with respect to E and A, under the pseudo-Newtonian potential 0i. The plot in the top left corner gives the 
variation for the regions of a single critical point only. The remaining three figures are for the multi-transonic region. Clockwise from the 
top right corner will give the variations for the inner critical point, the outer critical point and the middle critical point, respectively. 
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Figure 3. Variation of fl^ witli respect to A, for fixed values of £, in a polytropic flow. Eacli plot shows the birth of two critical points 
at a given value of A, which coalesce later at a higher value of A. 
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Figure 4. Variation of with respect to A, for fixed values of T, in an isothermal fiow. Unlike the case of polytropic flows, here the 
eigenvalues related to the inner critical point (a saddle point) continue to grow monotonically beyond the range of multi-transonicity. 



For the IPaczvnski fc Wiital ilQSOl) potential, 01 , the variation of Vi^ against £ and A is given in Fig. |5| The plot in the 
top left corner of this panel gives the dependence of Vi^ on E and A for the regions of the single critical point. The fiat portion 
is related to the outer critical point, while the hump is related to the inner critical point. Both of these points are saddle 
points. The vacant wedge-shaped area between these two regions is the region of multi-transonicity, which has been separately 
represented for the behaviour of the three distinct critical points in the remaining three plots in the panel. The top right 
corner shows the variation of the inner critical point, which is a saddle point. The bottom left plot gives the variation of the 
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eigenvalues for the centre-type middle critical point, while the properties of the outermost saddle-type critical point are given 
by the bottom right plot in Fig. |5| 

A common feature of all the three separate plots for the eigenvalues of the multi-transonic region is that they contain 
information on both accretion and wind. Here the lightly shaded regions (coloured red in the online version) in the plots 
indicate accretion. The criterion to distinguish accretion from wind has been that the flow rate through the inner critical 
point (a saddle point), A^in, has to be greater than the corresponding flow rate through the outer critical point (again a saddle 
point), TMout- The exact reversal of this requirement will give the wind region, which has been given by the darkly shaded 
portion (coloured blue in the online version) of each plot. These points can be better appreciated from the colour images on 
the archive version of this paper. 

Another property that one could discern from the surfaces in the {Q.'^,E,X) space is that the surface pertaining to the 
centre-type point has a minimum and then it rises, while the surface for the second saddle point has a miximum and then it 
dips. So the two of them taken together, they coalesce and "annihilate" each other at a large value of the angular momentum, 
and then all that is left behind is a single saddle point for the high angular momentum region of the parameter space. 

This contention is very illustratively borne out by the top left plot of Fig. |3 in which for all the three critical points 
has been plotted against A, for a suitably chosen value of f , in the case of a flow driven by the IPaczvhski fc Wiital lll980l) 
potential. Initially there is region of a single saddle point, followed by the birth of a centre-type point and another saddle point. 
Physically this is what it has to be, because with a centre-type point but without another saddle point, the flow solutions 
will all curl about the centre-type point, and there will be no means of connecting the event horizon with inflnity through a 
solution. To avoid this, and to make accretion a feasible proposition, the existence of the inner saddle point should be most 
crucial. The other three plots in Fig. |3 demonstrate this same feature for the other pseudo-Newtonian potentials. 

The whole situation is qualitatively si milar for isotherm al flows. The case for multi-transonicity in the parameter space 
of T and A has been discussed in detail bv lPas et alJ (1200,?^ . It is a simple extension from here to study the dependence of 
on T and A. This has been shown in Fig. |1] which gives the variation of 57^ against A, for fixed values of T. Once again 
it is evident that the multi-transonic region begins and ends, respectively, with the production and annihilation of a pair of 
critical points — one a saddle and another a centre-type point. A crucial difference for isothermal flows, however, is that 
beyond the multi-transonic region, for high values of angular momentum, the eigenvalues in relation to the lone saddle point 
increase monotonically in magnitude. 

A few points are worth emphasising upon once again at the end of this whole discussion, dwelling on both kinds of flows 
— polytropic and isothermal. Without encountering the difficulties of numerical integration, much predictive insight has been 
derived about the qualitative character of the flow, not the least of which, crucially, is the exclusive nature of the critical 
points, and in relation to that, important physical features of multi-transonicity itself. 



5 TIME-DEPENDENT STABILITY ANALYSIS OF STATIONARY SOLUTIONS 



Under the condition of hydrostatic equilibrium in the vertical direction, the time-dependent generalisation of the governing 
equations for an axisymmetric pseudo-Schwarzschild disc is given as 



dv dv IdP , A^ „ 

Ot or p or r'^ 

and 

9E 1 a ,^ , „ 

— -I--— E«r =0 
ot r or 



(20) 



(21) 



in which the surface density of the thin disc E, is to be expressed as S = pH JPrank et alJboO^) . Making use of Eq.@ and 
the polytropic relation P = Kp'' , Ea. l2H can be rendered as 



dt 



[P 



{7+l)/21 



+ 



7^3/2 0^ 



(7+l)/2 
p^ ' " V 



.3/2 



: 



(22) 



Defining a new variable / — p'^^'^^^^^vr'^^^ /^/Ij/, it is quite obvious from the form of Eg. 1221 that the stationary value 
of / will be a constant, /o, which c an be closely identified wit h the matter fiux rate. This follows a similar approach to 
spherically symmetric flows made bv lPetterson et all (|l98(l) and lThexms fc DavidI (|l992l) . For an axisymmetric disc, which 
has no de pendence on any angle variable, this approach has also been adopted for a flow driven simply by a Newtonian 
potential iRavll20o3l . The present treatment, of course, is of a more general nature, the disc flow being driven by a general 
pseudo-Newtonian potential, 4>{r). In this system, a perturbation prescription of the form v{r,t) = vo{r) + v'{r,t) and 
p{r,t) = po{r) + p'(r,t), will give, on linearising in the primed quantities. 



f' = fo 



p_ 

Po 



V 



(23) 
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with the subscript denoting stationary values in all cases. From Ea. l22L it then becomes possible to set down the density 
fluctuations p' , in terms of /' as 



2lL,n2 ( df_ 
dt fo [ dr 







with — 2(7 + 1) ^, as before. Combining Eas. (l231 and (1241 will then render the velocity fluctuations as 
dt fo \dt dr 

which, upon a further partial differentiation with respect to time, will give 

From Eg. 1201 the linearised fluctuating part could be extracted as 
dv' d 



d 






d 






dt 
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m] 



dt dr 



po 



(24) 



(25) 



(26) 



(27) 



with Cso being the speed of sound in the steady state. Differentiating Eg. 12711 partially with respect to t, and making use of 
Egs. 12411 . 1251 and 1261 to substitute for all the flrst and second-order derivatives of v' and p' , will deliver the result 

fo \ dt 



d_ 

dt 



vo df_ 
fo \ dt 



+ 



d 



dt 



df 
fo V dr 



+ 



d 



dr 



+ 



d 



dr 



^0(2 n2 2 \ df 

fo ^ ' dr 



= 



(28) 



all of whose terms can be ultimately rendered into a compact formulation that looks like 
9m {i^^d.f) = 

in which the Greek indices are made to run from to 1, with the identification that stands for t, and 1 stands for r. An 
inspection of the terms in the left hand side of Eq.JJHJ will then allow for constructing the symmetric matrix 

^ti-v _vo_ f 1 

fo \vo v'o - P^cio ^ 

Now the d'Alembertian for a scalar in curved space is given in terms of the metric g^^ by JVisseilll99d) 
1 



(29) 



(30) 



'd.4,) 



with g*"^ being the inverse of the matrix implied by g,,^ . Using the equivalence that i^" = g g*"^ , and therefore g - 
it is immediately possible to set down an effective metric for the propagation of an acoustic disturbance as 



?Gff 



Vo 



v'o 



(31) 
:det(f^''), 

(32) 



which can be shown to be entirely identi cal to the me tric of a wave equation for a scalar field in curved space-time, obtained 
through a somewhat different approach iVisseJ [19981) . The inverse effective metric, g^^, can be easily derived by inversion 
of the matrix given in Eg. 1321 . and this will give Vq — 0^c^o the horizon condition of an acoustic black hole for infiow 
solutions (^isscr 1998). 

A little readjustment of terms in Eg. 1281 will finally give an equation for the perturbation as 



9V n d 
— — + 2 — 

dt^ dr 



VO' 



df_ 

dt 



1 d_ 

dr 



vo 



( 2 
Vo [vo 



P Cso) 

' dr 







(33) 



whose form, as can be easily seen, remains invariant under any possible choice of a pseudo-potential. This is entirely to be 
expected since the potential, being independent of time, will have its explicit presence in only the stationary background fiow. 
Once again it is a simple exercise to check that an expression with the same form as Eg. 1)33^ could be derived for the case of 
isothermal fiows, but which will, however, have 13 = 1 and Cao as a constant. 

The perturbation is now fashioned to behave like a travelling wave, whose wavelength is constrained to remain small, i.e. 
it is to be smaller than any characteristic length s cale in the system. This kind of a treatment has been carried out before on 
spherically symmetric fiows JPetterson et al.lll980l) and on axisymmetric flows driven by a Newtonian potential llRavl2003l) . In 
both these cases the radius of the accretor was chosen as the characteristic length scale in question, and the wavelength of the 
perturbation was chosen to be much smaller than this length scale. In the present study, which is devoted to an axisymmetric 
flow driven on to a black hole by a pseudo-Newtonian potential, the radius of the event horizon could be a choice for such a 
length scale. As a result, the frequency, uj, of the waves should be large. A solution of the kind f'{r,t) = gui{r) exp(— ici;t) is 
used in Eg. 1331 . to give the expression 
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("0 



^2 2 N d ffc 



+ 



duo 1 d / 2 2 \ „. 1 dg^ /„. duo , 2\ ^ 

3uo— — tVoP Cso - 2lC<JDo -; 2lUJ— h (X) ffo; = 

dr vo cir ^ ' \ ar \ &r J 



(34) 



beyond which stage, mindful of the constraint that uj is large, the spatial part of the perturbation, g^ir), is first prescribed 
to be a power series in the form 



exp 



E 

n = -l 



(35) 



and this is then used in Ea. H,'i4^ . Following this, the three successive highest order terms in ld will be obtained as tj^, a; and {J* . 
Collecting all the coefficients of each of these terms, summing these coefficients up, and then separately setting each individual 
sum to zero, will yield, respectively, for u^, u and the conditions 
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and 



2 2 \ /dfc-lV 



21110- 
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dr dr 
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+ 2 



d/c_i d/ci 
dr dr 
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1 d 



d«o 

dr Ho dr 



(i;o/3^Cso 



dfco 
dr 



dfci 
2i«o-r^ = 
dr 



Out of these, the first two, i.e. Eas. H36|l and . will deliver the solutions 



and 



Wo ± /3CsO 



■ dr 



(36) 
(37) 

(38) 
(39) 



fco — — — In (/3uoCso) + constant 



(40) 



respectively. The two expressions above give the leading terms in the power series of g^(f). For self-consistency it will be 
necessary to show that successive terms follow the condition cij~"|/c„(r)| ^ cj"'""'"'^-' |fc„+i (r)j , i.e. the power series given by 
g^iir) converges rapidly as n increases. This requirement can be shown to be very much true, if the behaviour of the first three 
terms in fcn(r) is any indication. From Eqs. 1391 . (1401 and 1381 . respectively, these terms can be shown to go asymptotically 



as fc_ 



fco ~ Inr and fci 



given the condition that iio 



-5/2 



on large length scales, while Cso approaches its 



constant ambient value. Under the regime of high-frequency travelling waves, this will imply a;|fc_i| ^ |fco| S> a;~^|fci|, and, 
therefore, it should suffice here entirely to consider the two leading terms only, as given by Eqs. 1391 and 1401 . in the power 
series expansion of g^ir), as Eq. 1351 gives it. So with the help of these two terms, one may then set down an expression for 
the perturbation as 



/'(ni) 



y/3uoCso 



— iojt 

e exp 



«o ± /3cso 



■ dr 



(41) 



which should be seen as a linear superposition of two waves with arbitrary constants yl+ and Both these two waves move 
with velocity /3cso relative to the fluid, one against the bulk flow and the other along with it, while the bulk flow itself has 
velocity uq. With the help of Eg. 1241 it should then be easy to express the density fluctuations in terms of /' as 



,2 t'OPO 



/' 



/o V ^0 ± /3CsO 



and as a further step, with the help of Eg. 1231 . it should also be possible to set down the velocity fluctuations as 



±/3 



UqCsO 



/' 



with the positive and the negative signs indicating incoming and outgoing waves, respectively. 
In a unit volume of the fluid, the kinetic energy content is 



\ {po + p) {vo + v'Y 



(42) 



(43) 



(44) 



while the potential energy per unit volume of the fluid is the sum of the gravitational energy, the rotational energy and the 
internal energy, and is given by 



£pot = (po + p')(l)- (po + P') ^ + + '^'^ (^o"^) + (^o*^) 



(45) 
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where e is the internal energy per unit mass (iLandau fc Lifshit j|l987|) . In Eas. 1441 and 1451 . the first-order terms vanish on 
time-averaging. In that case the leading contribution to the total energy in the perturbation comes from the second-order 
terms, which are given by 

1 . ^2 



1 ,2 



■Sport = -^pov + vop + -^P (poe) (46) 

If the perturbation is constrained to be adiabatic, then the condition ds = can be imposed on the thermodynamic 
relation, de — Tds + [P/ p^^ dp, which will then give 



^2 

i£0 (47) 
Po 



Combining the results from Eas. H4H . 142^ . 1)43^ . 1)46^ and l|47|l . and upon time-averaging over the square of the fluctuations, 
will give 



2 1/32 



± "0 



(48) 



2 /(? («o±/3c,o)' 

with the factor 1/2 in Eg. 1481 arising from the time-averaging of the fluctuations. 

The total energy flux in the perturbation is obtained by multiplying fpert by the propagation velocity {vo ± /3cso) and 
then by integrating over the area of the cylindrical face of the accretion disc, which is 2-KrH. Here H is to be substituted from 
Eq.l^J, and together with the condition that H ^ r in the thin-disc approximation, and the fact that the speed of sound has 



a functional dependence on the density going as c^q = jKpQ ^ , an expression for the energy flux will be obtained as 



fo 



2/3 {M ± (3) 



(49) 



where M = M(r) — vq/cso, is the Mach number, which under a generic outer boundary condition, tends to vanishingly small 
values on large length scales of the accretion disc. All other factors in the expression of are constants. It is quite evident once 
again that the forms of both f pert and are invariant under the choice of any pseudo-potential, and therefore all arguments 
for the stability of the flow will be valid for any choice. Indeed, this should be especially true on large length scales of the 
flow, where all (j){r) will asymptotically tend to the Newtonian limit. 

As regards the question of the behaviour of the perturbation, it can be seen that for sub-critical solutions, the condition 
M{r) < /3 always holds everywhere for both the incoming and the outgoing modes of the perturbation, which, of course, 
implies that the total energy in the wave remains finite as the wave propagates. The perturbation manifests no unbridled 
divergence anywhere, and all sub-critical solutions are, therefore, stable. 

For the critical solution, it is evident that with the choice of the upper sign in Eg. 1)49^ . i.e. for incoming waves, which 
propagate along with the bulk flow, once again there is no growth behaviour for the perturbation to be seen anywhere and 
all stationary inflow solutions are, therefore, stable for this propagating mode of the perturbation. 

For the choice of the lower sign, i.e. for outwardly propagating radial waves, it is flrst of all very obvious from Eg. 1491 1. 
that there is no unrestrained growth of the energy flux of the perturbation in the region where M{r) < (5, and so the disc 
system remains stable in this sub-critical region. The case for stability is even better argued in the region of the flow where 
M{r) > 13. In this region, which has a finite spatial extent, all outwardly propagating waves will be swept away in finite time 
by the inward bulk flow right up to the event horizon of the black hole, and this physical e ffect will also ensure the stability 
of the flow in this region. This very same line of reasoning was adopted by ICarlicM (Il979ll to argue for the stability of the 
critical solution in the supersonic region of a spherically symmetric flow. As an interesting aside, it may also be mentioned 
here that this argument is well in consonance with the analogy of an acoustic black hole. Any disturbance within the sonic 
horizon for inflows, will not be able to propagate outwards through the horizon. Therefore, it can be established generally 
that in an inviscid pseudo-Schwarzschild system, all physically well-behaved inflow solutions will be stable under the influence 
of a linearised high-frequency radial perturbation, propagating either inwardly or outwardly. 



6 CONCLUDING REMARKS 

Saddle points are inherently unstable, and to make a solution pass through such a p oint, after startin g from an outer boundary 
condition, will entail an inflnitely precise fine-tuning of that boundary condition iRav fc BhattacharieQi2002il . Nevertheless, 
transonicity, i.e. generating a solution through a saddle point, is not a matter of doubt. The key to this paradox lies in 
considering a time-dependent evolution of the flow. In this paper, time-dependence in the flow has indeed been taken up, 
although in a perturbative manner, which has shown that all solutions of physical inter est are stable. So the selec tion of 
a solution through a saddle point will have to fall back on non-perturbative arguments jRav fc BhattacharieelEooi) . This 
is a relatively simple proposition for the case of a flow in the Newtonian framework. The difficulties against this line of 
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attack, however, are understandably much greater when the flow properties are studied in a general relativistic framework. 
Notwithstanding this fact, as a preliminary foray into this question, it should be important to have some understanding of 
the behaviour of the general relativistic flow solutions in their stationary phase portrait. This will be a matter of a subsequent 
investigation, for flows both on to a Schwarzschild black hole and a Kerr black hole. 
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